Direct experimental verification of quantum commutation relations for Pauli operators 
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We propose and demonstrate scheme for direct experimental testing of quantum commutation 
relations for Pauli operators. The implemented device is an advanced quantum processor that 
involves two programmable quantum gates. Depending on a state of two-qubit program register, 
we can test either commutation or anti-commutation relations. Very good agreement between 
theory and experiment is observed, indicating high-quality performance of the implemented quantum 
processor and reliable verification of commutation relations for Pauli operators. 
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Quantum theory associates each observable physical 
quantity with Hcrmitian operator acting on Hilbert space 
of states of a given physical system [T|. A fundamen- 
tal property of operators representing different quantities 
such as position, momentum or angular momentum of a 
particle is that they do not mutually commute, which 
gives rise to peculiar quantum effects like Heisenberg 
uncertainty relations. Although predictions of quantum 
theory have been corroborated by countless experiments, 
direct observation of the non-commutativity of the un- 
derlying operators has eluded us. Recently, testing of 
commutation rules for bosonic creation and annihilation 
operators based on the combination of single-photon ad- 
dition [2] and subtraction [3H5] has been reported [BHS]. 

Besides non-commutativity of bosonic creation and an- 
nihilation operators, commutation rules of Pauli opera- 
tors are also fundamental and important. The Pauli op- 
erators were originally introduced to describe Cartesian 
components of electron spin. More generally, they form, 
together with the identity operator /, a complete basis 
in the space of operators acting on a two-level quantum 
system, a qubit in the language of quantum information 
theory [5]. Pauli operators are both Hermitian and uni- 
tary, which can be represented by 2 x 2 matrices, 
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They satisfy fundamental commutation relations charac- 
teristic of Lie algebra su(2), 

[A, Y] = 2iZ, [Y, Z] = 2iX, [Z, X] = 2iY, (1) 

where [A, B] = AB — BA. Any two different Pauli oper- 
ators anti-commute, which means that 

{X,Y} = {Y,Z} = {Z,X} = 0, (2) 

where {A, B} = AB + BA. The nonzero commutator im- 
plies that, e.g., ZX ^ XZ, the overall operation depends 
on the order of X and Z. By combining Eqs. ([!]) and 
@ we find that ZX = -XZ = iY. It follows that, in 



contrast to the case of bosonic creation and annihilation 
operators [B], it is impossible to demonstrate the non- 
commutativity of Z and X by applying the two different 
sequences of operations ZX and XZ to some input single 
qubit state The output states ZX\tp) = iY\tp) and 
XZ\tp) — —iY\ip) differ only by an overall phase that is 
not directly observable. Even applying the operation to 
a part of an entangled two-qubit state does not help. 

In this letter, we report on the direct experimental ver- 
ification of commutation relations for Pauli operators. 
Our optical scheme combines two programmable quan- 
tum gates |10H12j and an auxiliary maximally entangled 
two-photon Bell state [131 [H] to directly implement the 
(anti)-commutator of two Pauli operators. We have com- 
pletely characterized the commutators by quantum pro- 
cess tomography [15T - H7] . Our work directly reveals the 
peculiar algebraic structure underlying quantum theory 
and represents realization of an advanced quantum infor- 
mation processor [10] - 

We work with optical qubits encoded into polarization 
states of single photons whose Hilbert space is spanned 
by the basis states \H) and \V), representing linearly 
horizontally and vertically polarized photon, respectively. 
In Dirac notation we have X = \H){V\ + \ V){H\, Y = 
i\V)(H\ - i\H)(V\ and Z = \H)(H\ - \V)(V\. As the 
Pauli operators are unitary, they can be deterministically 
implemented by (a sequence of) optical wave-plates |18j . 
Our main tool is a programmable quantum gate jTOl [11] 
where operation on signal qubit is controlled by the state 
of program qubit. The employed linear optical gate |121 
[THJ[20] consists of a polarizing beam splitter (PBS), where 
signal and program photons interfere, and the projection 
of the output program photon onto diagonally linearly 
polarized state \D) — -^(\H) + \V)). If the program 
photon is prepared in state \D), the polarization state of 
the signal photon will be unchanged and hence identity 
operation is applied. However, if we prepare program 
photon in orthogonal state \A) = ~^(\H) — \V)), the 
operation Z will be applied to signal photon. The success 
probability of the gate is 4 and does not depend on the 
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input state of signal photon. 

In our experiment we combine two programmable 
quantum gates with intermediate unitary operation U 
on the polarization state of signal photon (see inset in 
Fig. 1). By preparing the two program photons in the 
maximally entangled singlet polarization state \^^) = 
±(\A)\D) - \D)\A)) ee ^(\H)\V) - \V)\H)) we obtain 
the following transformation of state of signal photon, 
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which is, up to a constant prefactor, equal to the commu- 
tator Cz,u — [Z,U]. By varying U we can thus directly 
test various commutation relations. Moreover, by prepar- 
ing the two-photon program state in the triplet Bell state 
|$~) = ^(\D)\A) + \A)\D)) ee ±(\H)\H) \V)\V)) we 
realize anti-commutator of Z and U. 

The experimental setup is shown in Fig. 1. We first 
generate two pairs of entangled photons by spontaneous 
down conversion. The photons pass through the half- 
wave and quarter-wave plates (HWPs and QWPs) and 
are superposed on the PBSs (see Fig. 1) to implement 
the desired quantum gates. To achieve good spatial 
and temporal overlap, the photons are spectrally filtered 
(AAfwhw = 3.2 nm) and detected by the fiber-coupled 
single-photon detectors [13] . 

We have experimentally characterized the commutator 
Czjj by full quantum process tomography [T5HT7] for the 
following five different U: X, Y, ^±f, and H, 

where H = denotes a Hadamard operation. We 

have reconstructed the completely positive map x that 
describes the transformation of density matrix p, 

Pout = X(Pin) = Cz,uPmC zu - 

We can see that x contains complete information on the 
commutator C. According to Choi-Jamiolkowski iso- 
morphism [3TJ [52], x can be represented by a positive 
semidefinite operator on a Hilbert space of two qubits. 
This operator possesses an intuitive and appealing phys- 
ical interpretation: it is proportional to a density ma- 
trix of a two-qubit state obtained by sending through the 
channel x one part of a maximally entangled Bell state 
|$+) = ^(\HH) + \VV)). We have probed each commu- 
tator by six different input states forming three mutually 
unbiased bases {\H),\V}}, {\D),\A)}, and {|iZ),|£)}, 
where \R) and \L) denote right- and left-hand circu- 
larly polarized states. Each output polarization state was 
fully characterized by measurements in those three bases. 
From the collected data the operator x was reconstructed 
by means of a standard maximum likelihood estimation 
algorithm [2"5] . 

Examples of the results are shown in Fig. 2a where we 
plot the real and imaginary parts of the reconstructed 
X for three different U. For all five tested U the theory 
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FIG. 1: Experimental setup. Femtosecond laser pulses 
(394nm, 150fs, 76MHz) pass through two main BBO crys- 
tals (2mm) to produce two pairs of entangled photons with 
an average count of 2.6 x 10 4 /s. The mutual delays between 
path 2, 3, 4 are controlled by movable prisms Adi and Ad-2- 
We incorporate after each PBS a compensator (Comp.) to 
counter the additional phase shifts of the PBS. The photon in 
mode 2 serves as signal photon and the photon pair in modes 
3 and 4 represents the two-qubit program. The inset shows a 
conceptual diagram of the core parts of the scheme involving 
two programmable quantum gates (PQG) interspersed with 
another transformation U on the signal qubit. 



predicts that Cz,u = KV where V is a unitary oper- 
ation and K is normalization prefactor. Therefore the 
corresponding x should be proportional to the density 
matrix of a pure maximally entangled state. In par- 
ticular, for V — X we expect Xth = 8| X I' _ )(\1/ _ | and 
for U = Y we have xth = 8|* + )(* + |, where |* + ) = 
+ \VH)). The experimental results shown in 
Fig. 2a are in very good agreement with these theoretical 
predictions and the fidelity of the reconstructed commu- 
tators, defined as normalized overlap of x an d Xth, reads 
F x = 0.912 ± 0.008 and F Y = 0.873 ± 0.009. The opera- 
tor x is normalized such that Tr(x) = Tr(C*C) = 2\K\ 2 . 
However, the normalization factor K cannot be deter- 
mined solely from the tomographic data without some 
reference. 

We have therefore performed additional calibration 
measurements. We have introduced a temporal delay 
between the signal and program photons so that their 
wave-packets did not overlap on the PBSs and they be- 
haved as independent entities. To ensure the calibration 
data and signal data are obtained under identical cir- 
cumstances, we have carried out an independent calibra- 
tion for each U. With this calibration data at hand and 
taking into account that the success probability of each 
programmable quantum gate is \ we can normalize the 
data and fix Tr(%) and \K\. As can be seen in Table I, 
the experimentally determined factors K coincide with 
theoretical prediction within the statistical error. The 
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FIG. 2: Experimentally determined commutation and anti- 
commutation operations. Plotted are the reconstructed com- 
pletely positive maps \ characterizing the tested commuta- 
tion relations (a) and anti-commutation relations (b). For 
each case both real and imaginary part of the 4x4 Hermitian 
matrix \ i s shown. The matrices are normalized according to 
the experimentally evaluated K factors given in Tables I and 
II (see main text). 



fidelities of reconstructed \ are also shown in Table I. 
All fidelities exceed 0.8 which indicates good agreement 
between experimental observations and theory for all five 
tested commutators. 

The matrix \ actually characterizes the commutator 
only up to a phase factor. Indeed, two different commu- 
tators C\ and C*2 = e l ^C\ would be represented by the 
same x, because C\pC[ = C-2.pC\. However, the phase cfi 
does play an important role in the commutation relations. 
Although we can't directly measure the overall phase, we 
can verify the relative phase relations between two com- 
mutators. Consider the commutators Cz.x — 2iY and 
Cz.y = —2iX. From the measurements reported so far 
we can infer that, with high fidelity, Cz,x — e l ^ x 2Y and 
Cz,y = e l ^ Y 2X where the phases (fix, <fiy remain unde- 
termined. We next choose U = -^=(X + Y), By linearity, 
but without making any further assumptions, we have 
C ZtU = [Z,U] = \/2{e^Y + e^ Y X), hence C z ,u de- 
pends on (fix — 4>y- We have experimentally determined 
this commutator and the result is shown in the right col- 



umn of Fig. 2a. Notice the nonzero imaginary part of the 
matrix x- The reconstructed \ exhibits a high overlap (fi- 
delity 0.913 ± 0.009) with the maximally entangled state 
\^ v ) = -^(\HV) - i\VH)). One can easily check that 

= ^=(Y - A) (g) /|$+). The relative phase factor 
— 1 between Pauli operators A and Y is consistent with 
the theoretically expected relationship e l ^ x ~^ Y ^ 1 = — 1. 
The measurement thus corroborates the expected phase 
relationship between the commutators C'z,x and Cz,y- 

After testing the commutation relations, we have pro- 
ceeded to verify the anti-commutation properties of Pauli 
operators. For this purpose, we have changed the two- 
photon program state to \$~) and performed complete 
tomographic characterization of the anti-commutator 
A-z.u — {Z,U} for four different operators /, Z, H, 
and -^{Y ± Z). Examples of results of tomographic 
reconstruction are given in Fig. 2b. Since Azj = 2Z 
and Az.z = 21, the corresponding x are proportional to 
density matrices of Bell states |<I> - ) and \& + ), respec- 
tively. The experimental results shown in Fig. 2b are 
in good agreement with theory, as witnessed by high fi- 
delities of the reconstructed completely positive maps, 
F Az l = 0.910 ± 0.008 and F Az z = 0.897 ± 0.009. The 
reconstructed x representing Az,h and Az.z are very 
similar, which is not surprising because {Z,H} = \[21. 
The difference is only in the normalization of x that was 
determined from calibration measurements and reflects 
the different amplitudes of Az t z and Az.h- All four ex- 
perimentally evaluated normalization factors \K\ and fi- 
delities are summarized in Table II. 

Finally, we directly observed the anti-commutativity of 



TABLE I: The fidelities F and normalization factors K for 
the five tested commutation relations [Z, U] are listed. The 
normalization factors -ftf ca iib were determined from calibration 
measurements, K t h represent the theoretical predictions. 



u 


F 






X 


0.912 ±0.008 


1.98 ±0.03 


2.00 


Y+Z 

V2 


0.800 ±0.019 


1.40 ±0.04 


1.41 


Y 


0.873 ± 0.009 


1.98 ±0.03 


2.00 


X+Y 
s/2 


0.913 ±0.009 


1.95 ±0.03 


2.00 


H 


0.852 ±0.016 


1.39 ±0.04 


1.41 


TABLE II: 


The same as Table 


I but results 


for anti- 


commutation operations {Z,U} are 


shown. 




U 


F 


ifcalib 


Kth 


I 


0.910 ±0.008 


1.97 ±0.03 


2.00 


Z 


0.897 ±0.009 


1.95 ±0.03 


2.00 


Y+Z 
V2 


0.901 ±0.011 


1.36 ±0.03 


1.41 


H 


0.869 ±0.011 


1.39 ±0.03 


1.41 
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FIG. 3: Four-photon coincidence dips. We plot dependence 
of the four-photon coincidence counts on transformation U 
parameterized by the rotation angle of HWP. The dips di- 
rectly demonstrate the anti-commutativity of Z and X (a) 
and commutativity of Z operator with itself (b). The coin- 
cidence counts were measured for six input states \H), \V), 
\D), \A), \R), \L) and summed up. The squares are experi- 
mental data and the solid line represents best sinusoidal fit. 
The error bars are determined from Poissonian statistics. 



two different Pauli operators. Since e.g. {Z,X} = 0, we 
should not observe any four-photon coincidences when 
U = X and program photons are prepared in state |$ - )- 
By rotating the central half-wave plate by angle a we 
set U — cos(2a)Z + sin(2cv)X and measured the depen- 
dence of total number of coincidences on a, that in theory 
should be proportional to cos 2 (2cv). The results are plot- 
ted in Fig. 3a together with a fit by a sinusoidal function. 
We can clearly see the dip in coincidences at a w 45° 
which is a direct manifestation of the anti-commutativity 
of Z and X. The visibility of the dip obtained from the fit 
reads V = 84.6 ± 0.5%. We have also tested the relation 
{Z,Y} = 0. The results were very similar and are not 
presented here. Finally, we have directly checked that Z 
commutes with itself, [Z, Z] = 0. We prepared program 
photons in state l^ - ) and scanned a over the interval 
[-45°, 45°]. The result can be seen in Fig. 3b with the 
visibility V = 88.7 ±0.5%. 

The coincidence dips plotted in Fig. 3 are slightly 
shifted with respect to the theoretically expected posi- 
tions. We attribute this effect to the imperfections of the 
optical elements employed in the experiment. Another 
source of experimental errors is noise in the state gener- 
ation and the imperfect overlap of photons on the PBSs, 
which reduces the visibility of four-photon interference. 

In summary, we have devised and implemented a lin- 
ear optical scheme that enables direct observation and 
testing of quantum commutation relations for Pauli op- 
erators. In this way we can directly probe the non- 
commutativity of quantum operators corresponding to 
different physical quantities, which is one of the corner- 
stones of quantum physics. The demonstrated scheme 



also represents an advanced programmable quantum 
gate, where the type of operation (commutation or anti- 
commutation) is decided by the state of two-photon pro- 
gram register. By altering the program state, a whole 
class of operations can be realized, including all linear 
combinations of [Z,U] and {Z,U}. Thus, besides being 
of fundamental interest, our work may also find applica- 
tions in quantum information processing. 
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